Two channel orbital Kondo effect in quantum dot with SO(n) symmetry 
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A scenario for the formation of symmetry protected non-Fermi-liquid (NFL) Kondo effect (KE) 
with spin variable enumerating Kondo channels is suggested and worked out. In doubly occupied 
symmetric triple quantum dot within parallel geometry, the NFL low-energy regime arises provided 
the device possesses both source-drain and left-right parity. Kondo screening follows a multistage 
renormalization group mechanism: reduction of the energy scale is accompanied by the change of 
the relevant symmetry group from SO (8) to SO(5). At low energy, three phases compete: 1) under- 
screening spin triplet (conventional) KE; 2) spin singlet potential scattering; 3) NFL phase where 
the roles of spin and orbital degrees of freedom are swapped. 



I. MAIN RESULTS 

The physics of two-channel Kondo effect (KE) that is 
marked by non-fermi liquid (NFL) behavior at low tem- 
peratures has been with us for more than three decades^. 
Most attempts to realize it in complex quantum dots 
(CQD) failed because interference between two channels 
is virtually unavoidable, leading to channel anisotropy 
and to an instability of the two-channel fixed point. To 
remedy this instability, a suppression of interchannel co- 
tunneling is attempted, using special design of the CQD. 
Apparently, the most successful attempt is realized in 
double quantum dot (DQD)2, where the interference is 
suppressed by Coulomb blockade. Another design which 
allows one to (at least) approach the elusive two-channel 
fixed point was suggested in Ref^, based on a structure 
composed of triple quantum dot (TQD) in serial geom- 
etry (Fig. QJ,). Here, the strong Coulomb blockade in 
the central dot minimizes (but does not completely elim- 
inate) interchannel interference. 

Yet another approach for achieving NFL two-channel 
Kondo regime is to swap the roles of charge and spin 
variables, i.e. to treat orbital or charge fluctuations as 
pscudospin variables causing Kondo screening, whereas 
spin projection quantum numbers serve as different 
channels^—. The idea of accessing the NFL two-channel 
Kondo regime in "natural" two-level systems in heavy 
particles was subject to criticismZr— . In these systems, 
pscudospin-flip is a generic tunneling process with a char- 
acteristic (long) time itun, unlike real spin flip processes 
which are practically instantaneous. As a result the ul- 
traviolet cutoff for Kondo effect is the energy ~ H/tt un . 
This energy is of the order of the distance to the next 
excited level in two-well potential. The latter interval is 
much smaller than the energy scale Ef for "light" elec- 
trons. As a result the energy interval available for the for- 
mation of the logarithmic singularity is too narrow and 
the resulting Kondo temperature is very small Tk <C A 
where A is the depth of the occupied level in the well 
relative Ef- Thus the strong coupling regime remains in 
fact unattainable as far as two-level systems in heavy 



particles serve as pscudospin. 

Natural and artificial nano-objects provide their own 
mechanisms of two-channel Kondo tunneling assisted by 
pseudospin excitations. One such mechanism was pro- 
posed for a "quantum box" connected to a lead by a 
single-mode point contac t 10 ' 11 . In this case the relevant 
operator which logarithmically scales the crossover from 
the high temperature to the low-temperature region is 
the capacitance C(T) ~ hx(T/Tjc)- In another model 
with 5(7(3) dynamical symmetry the excited state with 
parity-degenerate rotational levels (m = ±1) may cross 
the m = level due to interaction with a bath and thus 
become the source of orbital KE with spin playing part 
of tunneling channel^. 

In this work we offer a relatively simple realization of 
over-screened orbital Kondo effect where the two-channel 
regime is realized by two spin projections. The proposed 
device is composed of a doubly occupied TQD in contact 
with two terminals within parallel geometry^. It will be 
shown that the present model is free of the shortcomings 
pointed out in Refs.— — and that the role of higher ex- 
cited levels is completely different. Experimentally, this 
configuration may be realized in triangular arrangement 
of vertical dotsH (see Fig [lb). The TQD is doubly 
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FIG. 1: Axially symmetric triple quantum dot in planar (left 
panel) and vertical (right panel) geometries. 
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occupied in the ground state, and only singly occupied 
states are involved in cotunneling processes. A state |A) 
of the isolated TQD (with corresponding energy E\), can 
be either \a)=\ ± |) or |A) = |00) (singlet), (triplet) 
for a TQD occupied by either N=l or N=2 electrons 
respectively (see below). 

The starting point is an Anderson-like tunneling 
Hamiltonian, H = Hb + Hd + -Stun where 



ff*=EEE ^4 fc ,c wfcCT , ff d = J2 ^a|a><ai 

i=(,r feo- a=e,o A=ct,A 

H tun = V2 ^(^Ac] efcCT |a')(A|+H.c.)+ff.c.(l) 



ikaa' A 

Here, assuming configurations having both left-right 
(Z-r) and source-drain (s-cQ symmetry, the operators 

c L CT =( c lfca ± cLfe^)/^ are even and odd combinations 
of source and drain electron operators. The tunneling 
amplitudes Va — (A\V\a' ,iekcr) between the dot and the 
leads depend on the two-electron state A of the TQDi&. 

The low-energy spectrum E\ of the isolated doubly oc- 
cupied TQD consists of two spin doublets Esi and two 
spin triplet Er% (i = g, u for even and odd combinations 
of I and r states), (see Appendix for detailed descrip- 
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FIG. 2: (Color online) Upper panel: Energy levels in dou- 
bly occupied TQD renormalized within a scaling procedure. 
Lower panel: Logarithmic energy scale for a crossover from 
SO (5) to two-channel KE as a function of scaling parameter 
r\ = In D /D. 

tion of the eigenstates and eigenvalues of the Hamilto- 
nian Hd). Within the energy scale of the bandwidth D 
(exceeding the width of this multiplet) the spectrum of 
the isolated TQD is characterized by SO (8) dynamical 
symmetry. The dynamical symmetry group is formed 
by linear combinations of the operators A AA = |A)(A'| 
generating the spectrum of the Hamiltonian Hd (seei^ 
for a regular description of dynamical symmetries) . Hal- 



dane renormalizationii— implies that all the levels flow 
downwards with scaling invariants, 



E* A = E A {D) - tt-^a HD/T A ). 



(2) 



Here Ta = ^PoI^aI 2 is the tunneling rate for the state 
| A), \n(D/Do) = r) is the scaling variable and Dq is the 
conduction bandwidth —Do/ 2 < Ek < Dq/2. The renor- 
malization rates T\/tt depend crucially on A, and the 
scaling trajectories E\(D) intersect as in Fig. [5] due to 
the inequality Tt 9 > Tg u > Tt„ > ^s g (see Appendix 
for details). Various scenarios of multistage Kondo ef- 
fect are possible, including that illustrated by the flow 
diagram of Fig. where the level E$ u crosses level Et u 
before the level Et "overtakes" both of them. There 
is a window of input parameters, where the orbital KE 
emerges due to nearly degenerate orbital doublet/spin 
singlet forming as a result of Haldane flow at D w D, 
where charge fluctuations are frozen. At this stage, an 
application of the Schrieffcr- Wolff (SW) transformation 
generates the effective spin Hamiltonian for N — 2 chan- 
nels, and further RG transformation follows Anderson's 
poor-man scaling procedure for renormalization of the 
exchange constants^. 

The full 5*0(8) spin Hamiltonian H^w is written iniS. 
From Fig. [2] we conclude that the poor-man scaling 
procedure should take into account the evolution of dy- 
namical symmetry along the chain 50(8) — s- 5*0(5) — s- 
orbital SU(2). To illustrate the key points of transfor- 
mation from spin KE to orbital KE, we consider the sim- 
plified picture, where only one of the two triplets is taken 
into account, i.e. write down the spin Hamiltonian per- 
taining to the part relevant to the KE for 50(5) multiplet 
composed of one triplet T u and two singlets S u and S g . 

H sw = Y,E Tu X T ^' T ^ + E S ,X S ^ S " + (3) 

H ri=u,g 
Jl^u ' &uu "F <^2-R*u ' S uu -|- ^(R^g ■ Sgu "F Rgii ' ^ug) • 

Here /i = 1,0, 1 arc projections of spin 5=1, the lev- 
els Ea are renormalized in accordance with $Z§. The 
group generators forming the 05 algebra are three vec- 
tors S M , R M , R = Riig + Rg 2 ] and the scalar A: 

5+ = V2{X lu0u + A°" Iu ), S* = X lulu — A IuIu , 
R+ = \/2(X luSu — A s " Iu ), i^ = -(A°" s " +A S "°"), 
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Sij are the components of local spin operators for 
even and odd partial waves of the band electrons, 

s y=5 Dfcfc' Ysaa' 4kaT~aa'Cgk>cr>, W are the components 
of the three Pauli matrices. 



The total 50(5) multiplet of width Ago (5) determines 
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the corresponding single-channel Kondo temperature 



Tki = D exp 



,16 



(5) 



provided D » T K i > A SO (s). Here ]\ = poJ l} po is the 
electron density of states on the Fermi level e F . This tem- 
perature is, however not universal: at the energy scale 
Aso(5) the fine structure of the multiplet deter- 
mines the Kondo scattering. 

A two-channel orbital KE (2CKE) is possible in the 
situation shown in Fig. [2j where at D < D the two singlet 
levels E u ,g, renormalizcd in accordance with Eq. ([2]), 
form an "orbital" doublet separated by a gap Ats from 
the triplet E Tu . The conditions for such level crossing are 
discussed in the Appendix. In this configuration, the SW- 
like Hamiltonian can be written in terms of pseudospin 
Pauli matrices T and f as, 

H mh = -K h % + A',, 7: ]Tr^ 

-Ki% ^{cl ka C e .k>a + cl M C oMa ) - K 2 I^T z . a . 



(0) 



kk' a fJ 

(see the Appendix for a definition of the parameters Ki 
Here 

% = , 



I = x SgSg + x SuSu , 

T a = C l,k<J C °> k '< 7 > T (J = C o,ko- Ce > k '< J > 
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Tz -,° ~ 9 ( ^2 C e,ka Ce ' k ' a C o,ka C o-k'a J 

\kk' kk' / 



(7) 



The first term in the Hamiltonian ^ is an analog of 
the Zeeman term in the conventional KE. Its origin is 
the avoided level crossing E$ u — E$ g = Kh arising in 
the course of Haldane renormalization due to weak in- 
terchannel hybridization in the leadsi^. The spin-flip 
processes are absent in the singlet states. As a result 
the spin degeneracy is symmetry protected and the spin 
index is converted into the channel number. Thereby 
we arrive at a desirable situation, where spin projection 
quantum number plays the role of channel index and the 
two (orbital) channels are identical. This is the orbital 
2CKE in an effective magnetic field The chan- 

nel isotropy is protected by spin-rotation symmetry of 
the singlet state. 

The second stage of the RG procedure, that starts at 
the energy ~ D ^> Ats includes inter-channel Kondo 
cotunncling ~ Ku , Kj_, and indirect virtual processes via 



excited triplet state |T„) "inherited" from the first stage 
of the RG procedure derived from the Hamiltonian . 
The full system of scaling equations encoding the orbital 
KE is presented in the Appendix. Here we present the 
reduced system including the parameters relevant for the 
scale T* characterizing the two-channel regime: 
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(8) 



Here the 3rd order terms on the RHS are retained in 
accordance with the general theory of two-channel KE^. 
The energy scale T* is enhanced due to contribution from 
the enhanced parameter j2- The coordinates of the cor- 
responding fixed point are 

1 + yr+ 2bi - ioni] 



K|| t = 
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(9) 



In this way we arrive at the two-channel Hamiltonian 
with pseudospin operator as a source of Kondo screening 
and spin indices enumerating screening channels. The 
Zeeman operator is relevant for the 2CKE, and its influ- 
ence on the scaling behavior in the nearest vicinity of the 
quantum critical point may be described within confor- 
mal field framework^. A peculiar feature of the orbital 
2CKE described here is the multistage renormalization 
of the parameters of the bare Anderson Hamiltonian fol- 
lowed by an appropriate modification of the dynamical 
symmetry shared by the spectrum of the doubly occupied 
TQD as displayed in (Fig. |2J}. 

The present scenario of orbital 2CKE implies pecu- 
liar behavior of various observables, such as the tem- 
perature dependence of the tunneling conductance G(T). 
It is essentially distinct from the analogous behavior of 
GIT) in the "conventional" spin 2CKE. First, in the 
weak coupling regime T 3> Ago (5^, Kondo co-tunneling 
in the orbital KE occurs according to the single chan- 
nel scenario. Two channel over-screening occurs only 
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at the strong coupling regime T <C Ago(5), where the 
crossover to NFL phase takes place. Second, in the or- 
bital 2CKE, the Kondo temperatures in the weak and 
strong coupling regimes are essentially different parame- 
ters, whereas in the spin 2CKE the two scales are nearly 
the same, T* = cxTk2, where a < 1 (see, e.g.,— ^). 
Therefore the two asymptotic regimes for G(T) are 



G(0)/G 



\n~ 2 (T/T K1 ), T»T K1 
1 + ^Jt/T*, T < T* 



(10) 



Numerical estimates of the RG parameters using 

G/G 




FIG. 3: (Color online) Temperature dependence of tunneling 
conductance (solid curve). Low-T orbital 2CKE and high- 
T 1CKE spin SO(5) asymptotics are shown by the dashed 
curves. The input parameters are specified in the text 

input values of the energy intervals At ,5^ =0.03, 
A Su ,t u =0.015, A TujSg =0.056 (all in units of£> ) show 
that the Haldane RG transformation stops at D/Dq = 
0.053 (the dependence -D(A) is neglected in this esti- 
mate). As is know n 15 ' 16 , the Kondo temperature for 
CQD displaying SO(n) symmetry is non-universal and 
depends on the width of the corresponding multiplet. 
In our estimate 7Vi(Ago(5))= 0.0054Do- On the other 
hand, the two-channel Kondo temperature Tk2 found 
from Eqs. ((8j) without 3-rd order terms on the RHS 
is Z/f2=0.018Z?o 3> Tki- The reasons for such strong 
enhancement are (i) the inequality Iifi(Ago(5)) <C 
Txi(0)m, (ii) the contribution of the enhanced param- 
eter j'2 to the first two equations in the system ijHJ: in 
accordance with our numerical calculations 32/ j% ~ 1.5, 
hence T* /T K1 > 1. 

This result may be compared with the situation aris- 
ing in doubly occupied DQD in the serial geometry with- 
out source - drain symmetry, where the difference be- 
tween the high-energy and low energy scales is related 
to the crossover from single channel KE to 2CKE due 
to asymmetry between the source/dot and drain/dot 
coupling^. The crossover mechanism 1CKE— > 2CKE 
is different in that case (two-stage Kondo screening of 
two spins S( S)( j) attached to two electrodes with essen- 
tially differing Kondo scales T Ks > T K Jl). The 1CKE 



arises at essentially lower temperature Tki <C Tk s due to 
overscreening of the remaining spin 1 /2 by the remaining 
electrons in the source and drain leads. Comparing the 
two mechanisms, we see that unlike^, our mechanism 
leads to enhancement of the Kondo scale in the crossover 
1CKE -> 2CKE. 

In conclusion, we demonstrate that the transforma- 
tion of the conventional high-temperature 1-channcl spin- 
Kondo effect into unconventional low-temperature 2- 
channel orbital Kondo effect, where spin enumerates the 
screening channels, is possible in doubly occupied par- 
allel TQD. The mechanism of this transformation is re- 
lated to the singlet-triplet level crossing (Fig. [2]). The 
necessary precondition for such crossover is the presence 
of at least two spin singlets in the low-energy spectrum 
of CQD, so that the minimal symmetry group is 50(5). 
In case of TQD with N = 2 the low-energy spectrum 
is formed by 50(8) two-triplets/two singlets multiplet. 
The theory is developed for the 50(5) case, but the gen- 
eralization for higher n is straightforward. Remarkable 
feature of this mechanism is the effect of enhancement of 
the 2CK Kondo temperature in comparison with that for 
1CK regime. This effect is robust because the former "in- 
herits" the Kondo cloud from the latter. The quantum 
critical point 2CKE-1CKE may also be achieved using 
the mutual disposition of the pairs -Et„. 9 and Es ug (Fig. 
[2]) governed by the Coulomb blockade energies (see Ap- 
pendix) as control parameters. With reasonable effort, 
the experimental techniques reported in Ref. (— ) can be 
modified and employed to test the present prediction. 
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Appendix A 

1. Eigenstates of doubly occupied symmetric TQD 

The Hamiltonian of the isolated isosceles TQD is: 

a—l,r,c er 
a i—l,r 

+WY J (dl a d ia +H.c). (Al) 



Here e a are the energy level positions in the central (c) 
and side (I, r) wells of TQD, Q a and Qi C are intradot 
and interdot Coulomb blockade parameters, respectively, 
W is the tunneling integral between the central and side 
wells. We will be interested in the completely symmetric 
case, ei = e r = e, Qi = Q r = Q < Q c and Qi c = Q rc 



- characterized by the scale Tk s , and the 2CKE regime The Hamiltonian (|A1I) can be diagonalized by using 
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the basis of two electron wave functions 
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V2 



Pi\exi r ) 
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\T U ) = a — (3i\ex lr ), 
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(A5) 



Pi = V2W/A, h = V2W/(A + Q), 



= V2W/(e c + Q C - Q lc - e), 



In this basis, the Hamiltonian (|A1|) is decomposed into 
singlet and triplet matrices, 
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2e + Q 
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2e c + Qc J 



(A3) 



where £ = £ r 



The low-energy multiplet of two-electron states is: 
2W 2 2W 2 AW 2 



E Sg = e 

Et u — £ 

Es n = e 
E T g = e, 
with A = e — f 



A A + Q £ C + Q 
2W 2 



A ' 
2W 2 
A + Q' 

-c Qic- 



(A4) 



= ^1-/3?, a u = ^1-/31. 



The cigenfunctions corresponding to the energies (|A4 



Thus the lowest state of the isolated TQD is Es g . It is 
seen from (|A4[) that the singlet and triplet states alter- 
nate: Es g < Et u < Es u < Et 9 ■ The effective RG proce- 
dure renormalizing the eigenvalues of quantum dot once 
it is attached the leads^i, has been generalized for multi- 
level TQD in Ref. O In accordance with this procedure, 
the levels E^ move downward as a function of the scal- 
ing parameter rj=lii(D / Dq) with different slopes oc T\. 
The tunneling rates Ta for even and odd states (|A5[) 
obey the following hierarchy: IV > > Tt u > ^S„- 
This means that multiple level crossing is possible in the 
course of RG evolution. Besides, the pairs Et u and 
E$ n arc subject to level repulsion provided the left and 
right tunneling channels are not completely independent. 

The flow evolution of each level stops at E\ w D. The 
value D is specific for each of the four levels. At this point 
the SW transformation leads to the Kondo Hamiltonian, 
and the Kondo-stage of RG procedure starts in accor- 
dance with the poor man scaling procedure^ The phase 
diagram is quite complicated, because the two singlet and 
two triplet levels can intersect in several ways depending 
on the values of the model parameters W, e, A, Q, Q c , Qi C . 
Among these scenarios are: (I) S=l Kondo regime cor- 
responding to the scenario where T g is the ground state, 
and the flow pattern of dynamical symmetries with the 
RG procedure is 50(8) -> 50(5) -> 50(4) -> 50(3) (cf. 
the case of TQD with N = 4 studied in Ref. M); (II) 
Absence of KE corresponding to the scenario where S g is 
the ground state; (III) Orbital KE with almost degener- 
ate S g , S u ground state. Here we focus on two competing 
phases, namely 50(5) configuration, which involves two 
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singlet states and one triplet state, and orbital SU(2) 
configuration, where two singlets form (quasi) degener- 
ate pair and the triplet state is involved as a relatively 
soft excitation above this doublet. 



The singlets S g and S u will be the lowest (renormal- 
ized) states in the SW limit, i.e., 

Es a = E Su - 2K h < E Ta , E Tu , when: 



connecting different states of the octet, 
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(A6) 
- Kh- 



Here K^/2 = M/ r is the indirect tunneling amplitude 
between the side dots via the central dot and the leads 
arising in the course of renormalizatiorJ^ [see also Eq. 
(|MT|) below]. 
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The spin operators for the electrons in the leads are 
introduced by the obvious relations 



s f - 2 E E c l kaT ' J ' J " 
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2. SO(5)-symmetry 



When the two singlets \S V ) and the triplet \T u ,/j) 
are almost degenerate, the exchange Hamiltonian can be 
written as a sum of spin and orbital parts: 



(A7) 



The first term in Eq. (|A7j) describes the singlet-triplet 
transitions: 



The Hamiltonian (|A8I) possesses an SO(5) symmetry. 
The group generators of the 05 algebra are the vectors 
S U ,R U from (|A10[) . the operators intermixing g- and it- 
states, namely, the vector R = Rug + Rgu, 



R + = V2(X luSg -X s " 1 "), 

R- = (R + V, 

R z = -(I°« s » + A 5 " "), 



(A12) 



and the scalar operators A interchanging g, u variables of 
the degenerate singlets: 



A = i(X SuS <> -X s <> Su ) 



(A13) 



H spin = Y^E Tu t u X t ^ +Y, E s, xSvS " 
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_l_ jT C . <, 1 J ST Tt .0 1 T ST (Tit 1 ) . a, +T?( 2 ).= "1 

* ° lu^u °u u lu ±x, u a u ^ u 2u V ug gu ' gu ug J 

Here Es g/u (D) = E$ T Aa and the effective exchange 
constants are 



T _ a 2 T V 2 
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fg T = ~^, (A9) 
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The spherical components of the vector operators S M , 
R u , Rug and R g 2 u are defined via Hubbard operators 



Additional vertices appear in the effective spin Hamil- 
tonian (|A8j) at the second stage of Haldanc- Anderson 
scaling procedure! 17 ' 18 As a result, the exchange part of 
the Hamiltonian ([3]) takes the form 

Hcot = Jl&u ■ S„ + J2R« • S u + J3(R^g > ■ s gu + Ilg 2 ) ■ s ug ) 
+ JiS U • Sg + J 5 R 

+ ^(R-lg ' s ug + R-2g ' s g«) + J?^u ' ( s gu + s «g) 

+ ^Rm ' Sg + Jl0R • Sg + JllRu • (Sg U + S U g) 

+ Ji2(Ri s •s s „ + R 2s -s ug ), (A14) 

where Ri; = X SlSr Ri, R 2i = R 2 X SrSl . The coupling 
constants in the Hamiltonian (|A14[) are subject to renor- 
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malization. Their values (IA9I) with 



J 1 (D) = Jl v , J 2 (D) = jf„ T , 



J 3 (7J) = J. 



ST 



Ji(D) = (i = 4-12XA15) 



are taken as initial conditions for solving the scaling equa- 
tions. These equations can be written in the following 
form: 
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d]5 
dhid 

dje 
dhid 

dj7 
dhid 

djs 
dhid 

dj s 
dhid 

djw 
dhid 

dj 11 
dhid 

dj 12 
dhid 



[fi + ft + if + 37 + ill + iii (k + iia) 



+ 



i3+i6+il + ii 2 2 



-[2(iii2 + jVin) + jVO'e + J12) - «W 5 ], 
-[33(71 + k) + 37(35 + iio) - ?7iir(i6 + in)], 
- [if + i? + ig + j'i 2 o + i?i + in(ie + 312) 

33 +36 +38 +312 

2 

-pilis + 37(33 + is) - m ;r .j 2 ], 

-[Mil +J4) - mi r j 3 ], 

(73 + is) (75 + iio) + (72 + ig)(i6 + j'12) 



5*7(71 +J4) + in (72 +ig)] ! 

-[38(71 + J4) + 57(75 + iio) - mirUu + 312)], 

-[2(j4j9 + jVj'll) + jV(i6 + J12) - "T-Zriio], 

-[25*45*10 + 37(33 + is) - "W 9 ], 
-[311(71 + J4) + 5*7(72 + ig)]i 
-[312(71 + i4) - ™/ris]- (A16) 
= Po^i ( ' = 1, 12), d = poD and m; r = 



where ji 

From equations (|A16j) . one deduces the Kondo tem- 
perature, 



'-K2 



D 1 



2\/2: 



1 



_) V2m^A17) 



ii + i2 + v / (.?i+i2) 2 + 2?1 
which transforms to Eq. (5) of the main text at mi r — > 0. 



3. Orbital KE 

The singlets S g and S u become the lowest renormalized 
states in the SW limit, i.e., E s =E S -2M ir < E T ,E T , 



when: 

2W 2 Pl ( 1 



2 W 2 
< -: ^ - Mu 



Pl+Pl \A ' e c + Q c -Q lc -eJ ' A + Q 
2W 2 Q 



< 



A(A + Q) 

2W 2 ((3 2 - (3 2 2 ) ( 1 



< 



\A e c + Q c -Q lc ~e 



(A18) 



In this case the two-channel orbital Kondo effect can 
be realized. The corresponding cotunncling Hamiltonian 
has the form: 

H olh = M lr (X uu -X™) 



a 2 V 2 



2 Sf — £ / — 

kk a 

al V 2 



E xSgSg iM c ^ 



k'a 



2 £ F - E t, 



/ , ^ ^o,ka °' k a 



(A19) 



fcfc'CT 



2 e_p 



with C( e)0 ) ife(T = (cj.jbff ± c r , fc(T ), and avoided crossing 
of the singlet states taken into account. The Hamilto- 
nian (|A19p can be rewritten in terms of pseudospin Pauli 
matrices / and r: 

H orb = -K h T z + K\\T Z E t z.t 

+^(r + ^+r-^ 

-KiT z ^2{c\ ka c e ^ a + cl ka c ^ a ) - K 2 I^2t z ^. 

kk' a a 

(see also Eqs. (6), (7) of the main text). The coupling 
constants are: 



(A20) 



K h = 2M lr7 



<, + < V 2 



K ± 



2 sf — e' 
agOtuV 2 
sf — £ 

a 2 - al V 2 



Ki = K, - 9 u 



4 e F 



(A21) 



The anisotropic Hamiltonian (|A20[) describes two- 
channel orbital Kondo effect in an effective "magnetic 
field" The coupling constant K^ remains unrenor- 

malized (because all the terms ~ t 2 contribute to Ki), 
but it affects the renormalization of the other coupling 
constants. 
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